svyset [pweight=wgt]
svy: ologit clim_conc i.lat_id##i.millz_id i.rep_id i.cath_id i.sex i.ed

margins lat_id,  predict(outcome(1))
margins lat_id,  predict(outcome(2))
margins lat_id,  predict(outcome(3))

margins lat_id, pwcompare(pv) predict(outcome(1))
margins lat_id, pwcompare(pv) predict(outcome(2))
margins lat_id, pwcompare(pv) predict(outcome(3))

margins millz_id, predict(outcome(1))
margins millz_id, predict(outcome(2))
margins millz_id, predict(outcome(3))

margins millz_id, pwcompare(pv) predict(outcome(1))
margins millz_id, pwcompare(pv) predict(outcome(2))
margins millz_id, pwcompare(pv) predict(outcome(3))

margins rep_id, predict(outcome(1))
margins rep_id, predict(outcome(2))
margins rep_id, predict(outcome(3))

margins rep_id, pwcompare(pv) predict(outcome(1))
margins rep_id, pwcompare(pv) predict(outcome(2))
margins rep_id, pwcompare(pv) predict(outcome(3))

margins cath_id, predict(outcome(1))
margins cath_id, predict(outcome(2))
margins cath_id, predict(outcome(3))

margins cath_id, pwcompare(pv) predict(outcome(1))
margins cath_id, pwcompare(pv) predict(outcome(2))
margins cath_id, pwcompare(pv) predict(outcome(3))

margins millz_id#lat_id, predict(outcome(1)) 
margins millz_id#lat_id, predict(outcome(2)) 
margins millz_id#lat_id, predict(outcome(3)) 

margins millz_id#lat_id, pwcompare(pv) predict(outcome(1)) 
margins millz_id#lat_id, pwcompare(pv) predict(outcome(2)) 
margins millz_id#lat_id, pwcompare(pv) predict(outcome(3)) 

